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ABSTRACT 



Previous studies have found that the width of gamma-ray burst (GRB) pulse 
is energy dependent and that it decreases as a power-law function with increasing 
photon energy. In this work we have investigated the relation between the energy 
dependence of pulse and the so-called Band spectrum by using a sample includ- 
ing 51 well-separated fast rise and exponential decay long-duration GRB pulses 
observed by BATSE (Burst and Transient Source Experiment on the Compton 
Gamma Ray Observatory). We first decompose these pulses into rise, and decay 
phases and find the rise widths, and the decay widths also behavior as a power- 
law function with photon energy. Then we investigate statistically the relations 
between the three power-law indices of the rise, decay and total width of pulse 
(denoted as 5 r , Sa and 5 W , respectively) and the three Band spectral parameters, 
high-energy index (a), low-energy index (/3) and peak energy (E p ). It is found 
that (l)a is strongly correlated with 5 W and 64 but seems uncorrected with 5 r ; 
(2)f3 is weakly correlated with the three power-law indices and (3)E P does not 
show evident correlations with the three power-law indices. We further investi- 
gate the origin of 5j — a and 5 W — a. We show that the curvature effect and the 
intrinsic Band spectrum could naturally lead to the energy dependence of GRB 
pulse width and also the 5d — a and 5 W — a correlations. Our results would hold 
so long as the shell emitting gamma rays has a curve surface and the intrinsic 
spectrum is a Band spectrum or broken power law. The strong <5 d — a correlation 
and inapparent correlations between S r and three Band spectral parameters also 
suggest that the rise and decay phases of GRB pulses have different origins. 

Subject headings: gamma rays: bursts — method: statistical 
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1. INTRODUCTION 

The origin of the gamma-ray bursts (GRBs) still remains unclear though it has been 
found for more than 40 years and many progresses have been made. Generally GRBs have 
complex temporal profiles with plenty of pulses and pulse is the basic element of light 
curves. The pulses for most bursts overlap each other, while there also exist well-separated 
pulses within some bursts. A lot of statistical studies about pulses have been made 
because they can provides us valuable clues to the radiation mechanism and underlying 
processes of GRB. In early statistical analysis, the width of GRB pulses was found to 
be energy dependent, i.e., the higher energies, the narrower widthes (e.g., Link, Epstein 
& Priedhorsky 1993). Using the average autocorrelation function to study the average 
pulse width, Fenimore et al. (1995) showed that the average pulse width of many bursts 
dependence on energy is well fitted by a power-law function with the power-law index about 
-0.4. The result was confirmed by later studies (e.g., Norris et al. 1996; Norris, Marani & 
Bonnell 2000; Nemiroff 2000; Crew et al. 2003; Norris et al. 2005). Peng et al. (2006) 
employed two samples consisting of 82 well-separated pulses to test the relation and found a 
power-law anti-correlation between the full pulse width and energy. A power-law correlation 
between the pulse width ratio and energy is also seen in the light curves of the majority of 
bursts in the two samples within the energy range of BATSE (Burst and Transient Source 
Experiment on the Compton Gamma Ray Observatory). Recently it is found that this 
power-law relation can be extended to X-ray bands (see, Zhang & Qin 2008, Zhang 2008) 
as well as X-ray flare (Chincarini et al. 2010). In addition, Zhang & Qin (2007) showed 
that the pulse peak time, the rise time scale, and the decay time scale on energy are also 
power-law functions. 

The origin of the dependence of the pulse width on photon energy is still unclear to 
date. It has been suggested that the power-law relation could be attributed to synchrotron 
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cooling (e.g., Kazanas, Titarchuk & Hua 1998, Chiang 1998; Dermer 1998; Wang et al. 
2000), which gives a slope of -1/2 between the pulse width and photon energy, similar to the 
observed -0.4. However, the synchrotron cooling is unlikely to be responsible for the relation 
due to its too short timescale relative to the pulse width. Under the assumption that the 
Doppler effect of the relativistically expanding fireball surface (or, in some papers, the 
curvature effect) is important, Qin et al. (2005) showed that, in most cases, this power-law 
relationship would exist in a certain energy range, and, within a similar range, a power-law 
relation of an opposite trend between the ratio of the rising width to the decaying width 
and energy would be expected for the same burst. Shen, Song & Li (2005) also considered 
the curvature effect and found the slope of the power-law relation to be -0.1~-0.2, less than 
the observed -0.4. These suggest that the curvature effect is an important factor to form 
the relation. 

The observed spectrum of GRB is usually a broken power-law form. The smooth 
connection of two power laws is the so-called Band function, which can described well the 
observed spectrum (Band et al. 1993). Three parameters in the Band function are used 
to characterize the spectral shape: the low-energy power-law index a, the high-energy 
power-law index j3, and the peak of the spectral energy distribution E p . It is found that 
the distributions of the power-law slopes of pulse width and photon energy obtained by 
Peng et al. (2006) and Zhang & Qin (2007) have large dispersions for different bursts, 
while the GRB spectra also vary dramatically for two different bursts. Whether the energy 
dependence of pulse is connected with the observed spectrum in some way or not? If it is 
true, then what does it imply for the mechanism of energy dependence of GRB pulse and 
whether does it supply some useful clues to the origin of the pulse rise and the decay phase? 
These motivate our investigations below. In Section 2, we present the sample description 
and pulse modeling. The analysis results are given in Section 3. We give the possible 
implications of the statistical correlations in Section 4. Conclusions and discussion are 
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presented in the last section. 

2. SAMPLE DESCRIPTION AND PULSE MODELING 

The sample we selected comes from Peng et al. (2010) observed by CGRO/BATSE 
with durations longer than 2 s, which contains 52 individual FRED pulses with the peak 
fluxes are greater than 1.0 photon cm~ 2 s" 1 on a 256-ms time-scale (for further information 
about the sample selection and spectral analysis, see Kocevski et al. 2003; Peng et al. 
2009a; Peng et al. 2010). 

In order to investigate the pulse temporal properties, we must model these pulses with 
a pulse model. A form proportional to the inverse of the product of two exponentials, one 
increasing and one decreasing with time, derived by Norris et al. (2005), is a good model 
to describe the GRB pulses. Thus we select the pulse model to model our selected pulses, 
which can be rewritten as 

I(t)=AXexp[-T 1 /(t-t a )-(t-t a )/T 2 ], (1) 

where t is the time since trigger, A is the pulse amplitude, t s is the pulse start time, T\ and 
r 2 are the characteristics of the pulse rise and the pulse decay, and A = exp[2(r 1 /r 2 )] 1 / 2 . 

Similar to Peng et al. (2006, 2009a) and Hakkila et al. (2008) we use the nonlinear 
least squares routine MPFIT and develop and apply an interactive IDL routine to fit all of 
the pulses, which allows the user to set and adjust the initial pulse parameters manually 
before allowing the fitting routine to converge on the best-fitting model via the reduced \ 2 
minimization. For each burst we require that the signal should be detectable in at least 
three channels (in this way, the relationship between the pulse width and energy can be 
studied). The fits are examined many times to ensure that they are indeed the best ones. 

We demonstrate the results with the fitting to GRB 950624b (BATSE trigger 3648), 
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in which three well-separated pulses were observed in four BATSE energy channels. The 
narrow distribution of the fitting x 2 values per degree of freedom (see, the right panel in 
Figure 1) indicates that the two-exponential model is sufficient to model the pulse light 
curves. 

The parameter values for all identified pulses are obtained, including the pulse peak 
intensity (A), pulse onset time (t s ), effective onset time (t e ff), and peak time (T pea k), as 
well as the two fundamental timescales (ti and r 2 ) (see Table 2 in Norris et al. 2005). 
The effective onset time, t e ff, is defined as the time when the pulse reaches 0.01 times of 
the peak intensity. Both onset times are relative to the burst trigger time. According to 
the fitting parameters we can obtain the measured pulse temporal properties including the 
pulse width w = Ary e = r 2 (l + 2 In A) 1 / 2 , the pulse rise width, r rise = \w(l — k), and the 
decay width r rise = \w{l + k). Fitting the pulse width, the pulse rise width, the pulse decay 
width (in the logarithm) per channel as a function of the geometric means of the lower and 
upper BATSE channel boundaries (using 300-1000 keV for channel 4) (Norris et al. 2005) 
shows the power-law relations indeed exist and the power-law indices are thus obtained 
(let 5 W , 5 r , and 5a denote the indices of the power-law relations between w, T rise , and Td eca y 
and energy, respectively). Therefore, the dependence of the three time quantities on energy 
can be parameterized by the power-law index. Example plots of the relations between the 
temporal properties and energy are illustrated in Figure 2 and the three power-law indices 
are presented in Table 1. 

In order to obtain the consistent time intervals between the light curves and spectra 
in the same pulse the time-integrated spectra are reanalyzed for each pulse of our sample 
based on the analysis of Peng et al. (2009a), which provided a detailed data description and 
spectral modeling of our sample. Only the so-called Band model (Band et al. 1993) is used 
to model the pulse spectra in this work. In the end there are 51 pulses that are included in 
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our analysis after removing one bad pulse spectrum with Band model in the certain time 
interval. The three spectral parameters for all of the pulses are also listed in Table 1. 
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Table 1. A List of the Burst Sample and Various Parameters. 



Trigger 8 W 5 r 5 d a E p (keV) 



563 


-0. 


694 


± 





.008 


-0 


.566 


± 





.042 


-0 


.729 


± 





.019 


-0.481 


± 


0.082 


-2 


.669 


± 





145 


184.500 ± 7.450 


907 


-0 


.542 


± 





.016 


-0 


.270 


± 





.049 


-0 


.647 


± 





.023 


-0.188 


± 


0.065 


-2 


.869 


± 





130 


188.400 ± 4.740 


914 


-0. 


.353 


± 





.039 


-0 


.408 


± 





.080 


-0 


.323 


± 





.020 


-0.848 


± 


0.191 


-2 


.475 


± 





.114 


100.900 ± 6.990 


973-1 


-0 


.176 


± 





.007 


-0 


.170 


± 





.009 


-0 


178 


± 





.006 


-1.102 


± 


0.024 


-2 


.084 


± 





.034 


298.500 ± 12.700 


973-2 


-0 


230 


± 





.019 


-0 


.176 


± 





.016 


-0 


.256 


± 





.024 


-1.431 


± 


0.044 


-2 


.264 


± 





.130 


254.500 ± 25.600 


999 


-0. 


.250 


± 





.028 


-0 


.108 


± 





.047 


-0 


.307 


± 


0. 


.023 


-0.971 


± 


0.052 


-2 


.010 


± 





.067 


388.800 ± 36.500 


1406 


-0 


379 


± 





.072 


-0 


.466 


± 





.082 


-0 


.352 


± 





.069 


-0.616 


± 


0.127 


-2. 


.138 


± 





.033 


121.700 ± 6.620 


1883 


-0 


.207 


± 





.062 


-0 


.165 


± 





.047 


-0 


.230 


± 





.070 


-1.294 


± 


0.039 


-3 


.772 


± 





.850 


291.000 ± 19.000 


1956 


-0 


163 


± 





.005 


-0 


.079 


± 





.008 


-0 


.215 


± 





.014 


-1.132 


± 


0.093 


-2. 


.385 


± 





.117 


144.300 ± 10.300 


2083 


-0 


.406 


± 





.010 


-0 


.278 


± 





.045 


-0 


.437 


± 





.024 


-1.326 


± 


0.039 


-3 


.533 


± 





191 


74.380 ± 1.180 


2138 


-0 


.408 


± 





.083 


-0 


.240 


± 





.010 


-0 


.474 


± 


0. 


111 


-0.315 


± 


0.151 


-2 


.933 


± 





.408 


222.900 ± 14.800 


2193 


-0 


582 


± 





.033 


-0 


.374 


± 





.025 


-0 


.571 


± 





.012 


0.384 


± 


0.064 


-2 


.898 


± 





.155 


274.000 ± 6.480 


2387 


-0 


.224 


± 





.009 


-0 


.143 


± 





.017 


-0 


.261 


± 





.020 


-0.372 


± 


0.048 


-2 


.482 


± 





.049 


158.400 ± 3.200 


2484 


-0 


.337 


± 





.010 


-0 


.199 


± 





.016 


-0 


.407 


± 





.022 


-0.310 


± 


0.123 


-3 


.287 


± 





.549 


178.100 ± 8.220 


2662 


-0 


.490 


± 





.019 


-0 


.370 


± 





.031 


-0 


.544 


± 


0. 


.016 


-0.700 


± 


0.129 


-2 


.838 


± 





.383 


178.800 ± 12.500 


2665 


-0. 


.799 


± 





.128 


-0 


.583 


± 





.075 


-0 


.890 


± 





.146 


0.292 


± 


0.271 


-2 


.856 


± 





.137 


87.720 ± 3.310 


2700 


-0. 


.223 


± 





Oil 


-0 


.534 


± 





006 


-0 


108 


± 





.012 


-1.361 


± 


0.044 


-3 


.005 


± 





.614 


249.800 ± 16.300 


2880 


-0. 


.426 


± 





.002 


-0 


.269 


± 





.082 


-0 


.462 


± 





.021 


-0.659 


± 


0.132 


-3 


.114 


± 





.406 


134.500 ± 6.410 


2919 


-0. 


.274 


± 





.041 


-0 


.222 


± 





180 


-0 


.296 


± 





.025 


-1.166 


± 


0.057 


-2. 


.253 


± 





146 


289.000 ± 28.500 


3003 


-0 


095 


± 





.014 


-0 


.058 


± 





.027 


-0 


.112 


± 





.009 


-1.122 


± 


0.034 


-2 


.053 


± 





.099 


489.600 ± 47.900 


3143 


-0 


.288 


± 





.009 


-0 


.147 


± 





.038 


-0 


.351 


± 





.030 


-0.916 


± 


0.500 


-2 


.032 


± 





.077 


101.700 ± 26.500 


3256 


-0 


.613 


± 





.003 


-0 


.653 


± 





.027 


-0 


.598 


± 





006 


-0.028 


± 


0.162 


-2. 


.939 


± 





.236 


145.800 ± 6.190 


3257 


-0 


578 


± 





.036 


-0 


.449 


± 





.030 


-0 


.613 


± 


0. 


.039 


-0.231 


± 


0.065 


-3 


.040 


± 





.226 


195.900 ± 5.160 


3415 


-0 


.193 


± 





.021 


-0 


.041 


± 





.013 


-0 


.254 


± 





.034 


-1.000 


± 


0.158 


-2 


.093 


± 





103 


182.500 ± 27.000 


3648-2 


-0 


528 


± 





.012 


-0 


.626 


± 





.060 


-0 


.479 


± 





.013 


-0.862 


± 


0.333 


-2 


.232 


± 





.064 


98.210 ± 10.000 


3648-3 


-0 


.528 


± 





.031 


-0 


.273 


± 





.027 


-0 


.625 


± 





.028 


-0.748 


± 


0.052 


-2 


.802 


± 





155 


217.200 ± 7.180 


3765 


-0 


231 


± 





.038 


-0 


.257 


± 





.165 


-0 


.227 


± 





.015 


-0.872 


± 


0.022 


-2 


.733 


± 





.101 


308.800 ± 7.750 


3875 


-0 


.425 


± 





.024 


-0 


.280 


± 





142 


-0 


.464 


± 





.063 


-1.245 


± 


0.333 


-2 


.974 


± 





.219 


56.380 ± 6.090 


3954 


-0. 


117 


± 





.013 


-0 


.253 


± 





.022 


-0 


.064 


± 





.023 


-1.153 


± 


0.059 


-1 


.927 


± 





.049 


295.500 ± 32.100 


4350 


-0 


.136 


± 





.052 


-0 


.389 


± 





159 


-0 


.047 


± 





.012 


-1.467 


± 


0.105 


-2 


.357 


± 





.420 


263.500 ± 62.500 


5478 


-0. 


.427 


± 





.006 


-0 


.225 


± 





.117 


-0 


.496 


± 





.033 


-0.287 


± 


0.107 


-2. 


.896 


± 





.265 


158.100 ± 6.410 


5517 


-0. 


.288 


± 





.070 


-0 


.230 


± 





.026 


-0 


.321 


± 





.094 


-1.156 


± 


0.154 


-2 


.667 


± 





.000 


149.700 ± 14.500 
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§ 3648 




Fig. 1.— Example plots of GRB 950624b pulse fits for BATSE channels (left panel) and the 
histogram of x 2 in our sample (right panel). 
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3. ANALYSIS RESULTS 

3.1. THE RELATIONS AMONG THE POWER-LAW INDICES 

Employing a sample consisting of 24 long-lag pulses Zhang & Qin (2007) studied 
the relations between the pulse temporal properties (width, rise width, decay width and 
peak time) and energy and found the pulse temporal properties are power-law functions 
of energy. In this section, we recheck the relations for the following two reasons. Firstly, 
the sample we used is much larger than that of Zhang & Qin (2007). Secondly, the FRED 
bursts are temporally and spectrally distinguished from long-lag bursts (Peng et al. 2010). 

Displayed in Figure 3 are the histograms of the three indices. The distribution 
parameters are listed in Table 2. From Figure 3 and Table 2 we find the distributions of 
these indices share the similar distribution width but they have large dispersions. The large 
dispersions imply that the energy dependence of the temporal properties may not be the 
same for different bursts. In addition, it is found that most of the power-law indices are 
much smaller than those of long-lag pulse derived by Zhang & Qin (2007). 

We also examine the relationships between the power-law indices of the temporal 
properties on energy. Figure 4 shows the relationships between them. The Spearman 
rank-order correlation analysis of the three quantities are listed in Table 3. The straight 
lines are fitted to the points: (1) 5 r = (0.11 ± 0.01) + (1.31±0.03) 5 W ; (2) 5 d = (- 0.07 ± 
0.06) + (1.01 ± 0.02) 5 W ; (3) 5 r = (0.09 ± 0.01) + (1.25±0.03) 5 d . It is found that the 5 d 
is strongly correlated with 5 W and the slope between them is close to 1. And so the two 
indices may be viewed as mutual surrogates. While the other two index pairs are obviously 
less correlated, which strongly indicates that the temporal properties (rise and the decay 
times) of the pulses do not evolve independently from each other, instead, their evolution is 
tightly coupled. 
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Table 1 — Continued 



Trigger & w 8 r 5 d a /3 E p (keV) 



5523 


-0 


269 


± 





.009 


-0.231 ± 0.018 


-0 


.290 


± 


0. 


.004 


-1.266 


± 


0.173 


-2 


.382 


± 





188 


129.200 ± 15.900 


5601 


-0 


.335 


± 





.024 


-0.246 ± 0.041 


-0 


.380 


± 


0. 


.030 


-0.725 


± 


0.057 


-2. 


.295 


± 





.076 


223.200 ± 10.300 


6159 


-0 


.155 


± 





.041 


0.035 ± 0.008 


-0 


.237 


± 





.052 


-0.907 


± 


0.446 


-2 


.260 


± 





.039 


73.200 ± 8.610 


6397 


-0. 


215 


± 





.008 


-0.116 


± 


0.014 


-0 


.259 


± 





.017 


-0.648 


± 


0.042 


-2. 


.572 


± 





.104 


192.500 ± 5.440 


6504 


-0. 


634 


± 





.037 


-0.475 


± 


0.005 


-0 


.687 


± 





.049 


0.067 


± 


0.099 


-2. 


.729 


± 





156 


157.200 ± 4.970 


6621 


-0 


112 


± 





.013 


-0.074 


± 


0.035 


-0 


.126 


± 





.009 


-1.086 


± 


0.074 


-2 


.424 


± 





.079 


128.900 ± 5.960 


6625 


-0 


.304 


± 





.001 


-0.253 


± 


0.009 


-0 


.331 


± 


0. 


.006 


-0.873 


± 


0.120 


-2 


.841 


± 





124 


76.900 ± 2.000 


6657 


-0 


.571 


± 





.051 


-1.017 


± 


0.258 


-0 


.521 


± 


0. 


.058 


-1.043 


± 


0.420 


-2 


.334 


± 





.159 


95.180 ± 13.100 


6930 


-0. 


.264 


± 





.016 


-0.171 


± 


0.021 


-0 


.314 


± 


0. 


.012 


-0.639 


± 


0.131 


-2. 


.294 


± 





.046 


95.180 ± 4.340 


7293 


-0 


.495 


± 





.014 


-0.281 


± 


0.062 


-0 


.585 


± 


0. 


.039 


0.313 


± 


0.089 


-2. 


.979 


± 





120 


161.300 ± 3.460 


7295 


-0. 


.410 


± 





.013 


-0.319 


± 


0.020 


-0 


.443 


± 


0. 


.022 


0.404 


± 


0.094 


-2 


.628 


± 





.180 


373.500 ± 16.500 


7475 


-0 


.109 


± 





.010 


-0.141 


± 


0.012 


-0 


,092 


± 





.010 


-1.315 


± 


0.066 


-1 


.979 


± 





.029 


149.900 ± 14.300 


7548 


-0 


227 


± 





.001 


-0.296 


± 


0.051 


-0 


.195 


± 





.026 


-0.836 


± 


0.104 


-2. 


.250 


± 





.000 


176.800 ± 11.100 


7588 


-0. 


.210 


± 





.024 


-0.423 


± 


0.026 


-0 


.134 


± 





.027 


-0.729 


± 


0.208 


-2. 


.710 


± 





.097 


71.290 ± 2.100 


7638 


-0 


486 


± 





.012 


-0.517 


± 


0.010 


-0 


.480 


± 





.013 


-1.298 


± 


0.389 


-2 


.576 


± 





.090 


54.170 ± 5.300 


7648 


-0 


.690 


± 





.197 


-0.826 


± 


0.292 


-0 


.647 


± 





.162 


-0.588 


± 


0.100 


-2. 


.356 


± 





.000 


201.100 ± 10.300 


7711 


-0 


272 


± 





.031 


-0.400 


± 


0.025 


-0 


.219 


± 





.032 


-1.241 


± 


0.051 


-3 


.124 


± 





.806 


211.200 ± 12.800 


8049-1 


-0. 


.208 


± 





.003 


-0.289 


± 


0.034 


-0 


.156 


± 


0. 


.018 


-0.391 


± 


0.104 


-3 


.826 


± 





.550 


164.300 ± 6.330 


8049_ 2 


-0 


324 


± 





.032 


-0.430 


± 


0.012 


-0 


.252 


± 





.063 


-1.408 


± 


0.072 


-3 


.096 


± 





.410 


157.400 ± 12.100 



Table 2. A List of the Distribution Parameters of the Three Power-law Indices. 



Power-law indices Mean Median a (modeled with a Gaussian profile) 



S r 



-0.35 
-0.31 
-0.37 



-0.31 
-0.27 
-0.32 



-0.32 ± 0.17 
-0.27 ± 0.17 
-0.33 ± 0.23 
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Fig. 2. — Example plots of the relations between the observed pulse temporal properties and 
energy for our selected samples, where the open circle, the diamond, the triangle and the 
square represent the pulse width, the pulse rise width, the pulse decay width, respectively. 
The power-law relations between them are evident. 

Table 3. Correlations of the power-law index pairs. 



Power-law indict 



Rs 



Ps 



S w - 5 d 

<5r- S d 



0.69 2.63 x 1CT 8 
0.96 2.93 x 10~ 28 
0.49 2.84 x 10~ 4 



Note. — Rs and Ps denote the Spearman 
rank-order correlation coefficient and signifi- 
cance, respectively. 
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Fig. 3. — Histograms of the power-law indices 5 W (a), 8 r (b), and 5d (c) obtained by fitting the 
pulse width, rise width, and decay width and energy with power-law functions, respectively. 
The dashed lines are the best fits by the Gaussian functions. 
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Fig. 4. — Relationships of the three power-law indices 5 r vs. S w (a), 5a vs. 5 W (b), 5 r vs. 5d 
(c). The dashed lines are the regression lines. 



Table 4. Correlations of the Three Power-law Indices versus Spectral Parameters 



Parameter 


Rs 


Ps 


Parameter 


Rs 


Ps 


Parameter 


Rs 


Ps 


8 w -a 


-0.57 


1.41 x 10~ 5 




0.39 


4.42 x 10" 3 




0.20 


1.05 x lO" 1 


S r -a 


-0.24 


8.99 x 10~ 2 


5 T -P 


0.34 


1.37 x 10" 2 


6 r -Ep 


0.19 


1.77 x lO" 1 




-0.62 


1.04 x 10~ 6 


S d -I3 


0.32 


2.22 x 10" 2 


8 d -E p 


0.21 


1.44 x 10" 1 



Note. - 



Rs and Ps denote the Spearman rank-order correlation coefficient and significance, respectively. 
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3.2. THE RELATIONS AMONG THE POWER-LAW INDICES AND THE 

BAND SPECTRAL PARAMETERS 

Since the power-law indices reflect the spectral evolution of GRB pulse we wonder 
whether they are related to the spectral parameters or not. With the temporal and spectral 
parameters obtained in the previous section, we can examine the relationships between 
the power-law indices and the spectral shape parameters, low-energy index a, high-energy 
index j3, and peak energy E p . Illustrated in Figure 5 are the scatter plots of the power-law 
indices versus a. The Spearman rank-order correlation parameters are listed in Table 4. 
The regression analysis gives the best-fitting lines: (1) 5 W = (—0.63 ± 0.01) + (0.40 ± 0.01)a; 
(2) 5 r = (-0.70 ± 0.02) + (0.50 ± 0.02)a; (3) 5 d = (-0.62 ± 0.01) + (0.39 ± 0.01)a. 

The strong anticorrelated relations are identified between 5d and a as well as 5 W and 
a. However, there seem no correlation between 5 r and a and the dispersion is much larger 
than the other two parameter pairs. In addition, the slope and intercept of 5 r and a has 
apparent difference compared with that of 5d and a. A possible interpretation of the 
phenomenon is that the mechanism causing the dependence of the rise width on energy 
might be different from the other temporal properties on energy. Another possible reason is 
that the mechanism of producing the rise phase is different from that of the decay phase. 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

a a a 



Fig. 5. — Power-law indices 5 W (a), 5 r (b) and 5d (c) vs. low-energy power-law index a. The 
dashed lines are the best fitting lines. 
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Fig. 6. — Power-law indices 5 W (a), 5 r (b), and 5d (c) vs. high-energy power-law index 
The dashed lines are the best fitting lines. 
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Fig. 7. — Power-law indices S w (a), 5 r (b), and 5d (c) vs. peak energy E p . No correlation is 
apparent. 
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Examining of the power-law indices versus the other two shaped parameters (3 and E p , 
shows that all of the power-law indices are weakly correlated with the (3 (see, Figure 6 and 
Table 4). However, no apparent correlations are found between the three power-law indices 
and E p (see, Figure 7 and Table 4). 

4. IMPLICATIONS OF THE STATISTICAL CORRELATIONS 

The origin of the found correlations between the power-law indices and the Band- 
spectrum parameters could be related to the GRB pulse formation mechanism. The 
mechanisms of the rise phase and decay phase of the pulse are different. The rise phase of 
the pulse is believed to be related to gamma-ray active time of emission region, while the 
decay phase originates from the curvature effect. The found strong 5 a — a correlation and 
weak S r — a correlation also suggest the different origins of the rise phase and decay phase. 
It is interesting to speculate the physical origin of the strong 5^ — a correlation. Since 5^ is 
a quantity describing the decay phase, we consider the effect of the curvature effect on the 
energy dependence of the pulse. 

4.1. The energy dependent of pulse width due to the curvature effect 

In order to study the energy dependent of pulse width due to the curvature effect, we 
first calculate the decay profile of a pulse from a thin spherical surface expanding with 
relativistic speed. Define an emissivity 

fE> = Y!j(t'-t' )6(R>-R> ), (2) 

where YI'e' * s ^ ne surface brightness per unit photon energy interval in the comoving frame. 
Since the observed GRB spectrum is broken power law form, it is not unreasonable to 
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assume YUe' t° t> e the following form: 



y' _ y' J ( e '/ e p) 

~ ^° I {E'/E p) 



(3) 



{E'/E p ) 1+a E'<E' p 
) 1+ ? E>>E' p 

where V^ is the normalized coefficient, E' p is the intrinsic cutoff energy (comoving frame), 
and a and j3 are the low and high energy photon number spectrum slopes, respectively. As 



we know, more generally GRB spectrum is described by the so-called Band function: 



( 



E' 



lOOfceV 



(a-P)E' 



\ 1+a 

) exp 



E'(2+a) 



E' < E' p (a-p)/(2 + a) 



^P(/3-«)(lo£y) 1+/3 E>>E>(a-P)/(2 + a) 



(4) 



(2+a)100fceV 

However we will find below, the broken power law spectra gives more explicit physical 
meaning than the Band function. 

The observed photon flux density can be given by 

i r . i 



F F = 



4nD 2 



J E dV 



2D\ 



i' E ,D 2 R 2 dRdn, 



(5) 



Here D = l/[r(l — j3rfj)] is the Doppler factor and the transformations ]e = D 2 j' E , and 
E = DE' are used. The photon flux per unit photon energy can be given by 

F E 



N e (T) 



E 



(6) 



The observed time is T = t — Rfi/c and T = was chosen as the time of arrival 
at the observer of a photon emitted at /i = 1, t — t and R = R . Thus we have 
to = R /c and D \ R=Ro = 2T/[1 + (T/T ang )\ where T ang = R /(2T 2 c). Note (3 r = 1 - 1/2T 2 
is used in the above derivation. Using dfi = cdt/R, di! = dt/T, dR! = TdR, and 
6{t - t' )8(R' - R' ) = 5{t - t )5{R - R ), one can get 



N E (T) 



Fe 
E 

= N D 2 



{E'/E' p ) a+1 E'<E' p 
(E'/E' p Y +1 E>>E' p 



(7) 
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= N D 2 




(8) 



(^r 1 E>DE> 



where N = Sq_R c/ \2D 2 l E) is the normalized coefficient. One can find, from the above 
equation, the observed photon flux behaves as N E oc (1+T /T ang ) a ~ l or N E oc (1+T /T ang Y^ 1 
for a given observed photon energy E. The observed photon flux within an energy interval 
of Ei — E 2 can be written as 



The decay phase of light curve in GRB is usually believed to be the result of curvature 
effect, i.e., the photons at higher latitude arrive at the observer more later with lower flux 
due to the spherical emitting area and Doppler effect. If we consider the light curve for a 
given observed energy, then the comoving (or intrinsic) photon energy contributing to the 
given observed photon energy moves toward higher energy end with the decay of light curve 
due to the Doppler effect. The movement would experience a break so long as the given 
energy is less than E p since the intrinsic spectrum is a broken power law form. The light 
curve will first decay with N E oc (1 + T/T an9 ) a_1 , and then with N E oc (1 + T/T ang Y~ l 
so long as the given observed energy satisfies E < E p . The transition time, defined as T ep , 
is the time when the intrinsic photon energy contributing to the flux at E just moves to 
E' p . However, usually the full width half maximum (FWHM, defined as T1/2) of a pulse or 
the width when the flux decays to 1/e of peak flux (defined as Ti/ e ), as used in this paper, 
are considered to be the widthes of the pulse. These artificial definitions of pulse width 
could lead to the result that the energy dependence of pulse width with energy does not be 
measured though it can present. The reason is as follows. Following the the definition of 
T1/2, it can be given by Ti/ 2 = [2 _1// ( Q_1 ) — l]T ang . Consider an observed photons with energy 
E < E p . The intrinsic contribution to it is initially from photons of E' a = E a /(2T) < E' p . 
When the photon energy increase from E' a to E' p as the decay of Doppler factor in high 



N(T) 




(9) 
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latitude region, the Doppler factor is D^ a = E a /E' p = 2Y / (1 + T ep /T ang ) and the time is T ep , 
so T ep is solved to be T ep = (2TE' p /E a - l)T ang . If T ep > T 1/2 , i.e., a < 1 - ln2/ln(E p /E a ), 
then the FWHM is only determined by the low energy section, otherwise, the FWHM 
of light curve will be affected by the high energy part of the intrinsic spectrum. If two 
photons, say E a and Et>, both satisfy a < 1 — ln2/ln[E p / E a (b)}, then we can't measured 
the energy dependence of pulse width. Only if one of the photon energies or both of them 
satisfy the condition a > 1 — ln2 / ln[E p / E a (b)}, the energy dependence of pulse width would 
be measured. If we consider the time Tu e , then T ep < Ty e gives 

a > 1 - l/ln(E p /E a ). (10) 

In this condition, the closer to E p the photon energy (i.e., the higher energy), the narrower 
the pulse width of this energy, shown in Figure 8. This naturally explains the energy 
dependence of pulse width. If two energy bands are both more energetic than E p , then the 
pulse widthes in these two bands will be the same, i.e., no pulse broadening with energy 
presents in current consideration. 

BATSE bursts, the typical E p is ~ 200 keV, the typical low energy slope is —1 and four 
energy channels are (25-50, 50-100, 100-300, and 300-1000 keV) (Preece et al. 2000). The 
condition (|T0|) give E p > 121 keV for the typical parameters. Therefor for energy channels 
1 and 2, the condition is not satisfied and thus the pulse broadening in low energy band will 
not be measured in the two channels. However, it is so only for typical parameters. For a 
specific burst, a > 1 — l/ln(E p /E a ) may be satisfied. Channel 3 crosses the typical E p , while 
channel 4 is larger than it and thus only decays with Ne oc (1 + T /T ang Y~ l . Thus the pulse 
width will generally decrease from energy channel 4 to channel 2. In current consideration, 
our results indicate no broadening measured in channels 1 and 2 for the bursts with typical 
parameters. However, there are two possibilities to lead to the broadening in channels 1 
and 2. 
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One possibility is that the shell of internal shock is not a spherical symmetry, which is 
proposed by some authors, e.g., Fenimore & Sumner (1997) and Kocevski, Ryde & Liang 
(2003). Kocevski, Ryde & Liang 2003 found that only ~ 40% bursts are consistent with 
the spherical curvature. Their results point to a picture that the shell is a prolate geometry 
for a large number of bursts. In this case, the Doppler factor is decreasing faster than 
the spherical geometry so that the condition of (ITUj) is easier to satisfy. Thus the prolate 
shell could lead to the broadening in channels 1 and 2. Another one is that the intrinsic 
spectrum can be evolving. If the intrinsic spectrum is evolving, e.g., E p decreasing with 
time, then the condition (fTUj) can be satisfied so that the pulses in channels 1 and 2 also 
broaden with energy. However, we do not know how the intrinsic spectrum evolves and 
how fast it does. The curvature effect naturally gives rise to an observed spectral evolution 
due to the Doppler effect, which gives an E p decay with T~ x for a spherical curvature. If 
intrinsic spectral evolution gives a steeper observed one than T , then it will be masked 
and thus yield to T~ l given by curvature effect. 

The above analysis assumes the intrinsic spectrum is broken power law, so the light 
curve broken of decay phase at T ep can be seen clearly (see Fig. [8]). Actually the GRB 
spectrum is generally smoothly connected at the broken energy, namely, the Band spectrum, 
and the observed light curve is usually from an energy interval, not a single energy, so the 
observed decay phase is usually smooth. Note that total pulse width in GRB is dominated 
by decay phase. Thus the curvature effect+Band spectrum we consider here can provide a 
natural explanation to the energy dependence of GRB total pulse width, not only the decay 
width. 
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Fig. 8. — The light curves for different photon energy. Arrows mark the positions of T ep . 
The dashed lines label the positions of half and l/e maximum flux. Line a decays as (1 + 
T/Tangf- 1 , lines b, c, d and e first as (1 + T/T ang ) a - 1 and then as (1 + T/T ang f" 1 and line f 
as (1 + T jT ang ) a ~ x . Lines b and c show energy dependence of width of decay while the lines 
d and e do not due Ti/ 2 < T ep if considering FWHM as the decay width. 
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4.2. 5d — a relation 

We have shown the Band spectrum+curvature effect is a reason of pulse broadening 
with photon energy. Different Band function parameters will cause the differences of the 
degree of pulse broadening. However, the effects of the three Band function parameters, 
a, (3 and E p , on the degree of pulse broadening are different. The main effect seems to 
come from the low-energy slope. The reason is the low-energy section is usually flatter than 
the high one, i.e., a > {3. Thus the decay, N E oc (1 + T/T ans ) a_1 , is also more shallow 
than Ne oc (1 + T/T an£ ,)^ _1 with the pulse decay width dominated by a, unless the given 
energy E > E p . One can come to the conclusion that the larger a, the more significant 
the broadening with energy. Suppose two extreme cases: a = (3, i.e., the spectrum is a 
single power law; a = 1 and (3 = oo, i.e., the spectrum is composed of a horizon line and 
the vertical line. The pulse width would be the same for different energy bands in the 
former case, while the broadening of pulse with energy is the most significant in the latter 
case. Other cases lies between the two cases. It is difficult to find directly whether this 
is consistent with the observed anticorrelation in the earlier sections. We thus make a 
calculation to test that. 

In order to compare with the observations, we calculate the count rate within four 
energy channels (25-50, 50-100, 100-300, and 300-1000 keV) with equation EJ Further 
we calculate the decay width of pulse within the four energy channels, fit linearly the 
width-energy relation and then find the power law slope, 8^. Figure 9 shows our theoretically 
calculated a — 5<i relation. Obviously, a and 5d are roughly linearly anti-correlated, which 
is consistent with the statistical results of the observations. The fitted slope in Figure 9 is 
0.3, also roughly consistent with the statistical results. This suggests that the curvature 
effect -l-Band spectrum could lead to the observed 5d — a correlation. Strictly speaking, 
5^ and a are not linearly correlated for the theoretical results (still see Figure 9), not 
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completely consistent with the observations. However, in our calculations, some plausible 
values of the parameters, for instance the Lorentz factor and the peak energy, are used, 
disregarding the variety of the parameter value for different bursts, which could lead to the 
inconsistency Further whether the relation in physics is linearly anti-correlated or other 
forms is not important. The observed relation results from the sum of intrinsic physical 
factors, such as the Band spectrum, curvature effect and some artificial factors , such as 
the energy channel selection of BATSE and the pulse width definition, which is difficult 
to figure out one by one. Here we only focus on the Band spectrum and curvature effect. 
What we concern here is the effect of the combination of them on the pulse. We find that 
the curvature effect+Band spectrum can indeed result in such an anti-correlated trend of 
a - 5 d . 

5. Discussion and Conclusions 

In this paper, we have selected 51 well-separated long-duration FRED GRB pulses from 
the available sample of BATSE and analyzed the pulse temporal and spectral parameters 
with the aim to search for their connections. Our analysis first shows that the pulse width, 
pulse rise width and pulse decay width with energy all scale as a power law with the 
power-law slopes of 5 W , 5 r and 5^, respectively. In addition, the slopes are correlated to each 
other but the correlation between 5 W and Sd is much stronger than that of 8 W and 5 r as well 
as 8 r and 5d, which may be due to the pulse width is dominated by the decay phase and the 
rise and decay phase have different origins. 

Then we investigate the relations between the power-law indices and the three Band 
spectral parameters. Our results show that the power-law indices, 5 W , 5d are strongly 
correlated with the low-energy index a and much less correlated with high energy slope j3. 
The 8 r weakly correlates with (3 but does not show apparent correlated with a. 
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Fig. 9. — Sd — a relation. The solid lines are for 4 channels of BATSE considered, while the 
dashed lines for 1, 2 and 3 channels case. The latter is considered due to the fact that the 
fitted 5 r , 5d and 5 W in statistical part of the paper are obtained mainly using channels 1, 2 
and 3. 
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We do not find that there are apparent correlated relations between the peak energy 
E p and three power-law indices based on our current sample (see, Figure 7 and Table 4). 
However, we think the demonstrations that 8 W , 5 r , and 5d are independent of E p may be 
weak only based on our current sample. Several reasons may cause the case. Firstly, it is 
well known that the BATSE E p distribution is narrow (e.g. Band et al. 1993, Mallozzi 
et al. 1995, Schafer 2003), and thus it can be hard to clearly identify the correlations in 
data drawn from a narrow distribution. Secondly, the sample is likely biased toward low 
-Ep-valued pulses. The sampled pulses have been selected on the basis of long durations and 
bright peak fluxes. Because of pulse property correlations, such pulses tend to be softer than 
the mean (Hakkila and Preece, 2011). This pulse selection bias can have contributed to the 
narrowness of the sample. Lastly, some of the highest E p pulses in this sample may have 
been incorrectly identified as single pulses, when in actuality they represent merged pulses. 
If the hard-to-soft pulse evolution demonstrated by Peng et al. (2009a) is normal, then a 
pulse made from two merged pulsed will experience a re-hardening when the second pulse 
kicks in (Hakkila and Preece 2011; Ukwatta 2011), which can give these pulses higher E p 
values than the other single pulses to which they are being compared. Due to the bad fits 
to those merged pulses with given pulse model we are not sure how the merging of pulses 
affects the 5 values from the data analysis. But we can give some reasonable deductions. 
Firstly, it seems that the smallest 5 value from a merged pulse occurs when pulses overlap 
completely (i.e. when they peak at the same time). Otherwise, the pulses should appear to 
be longer in all energy channels, and the pulse separation will make a larger contribution to 
the pulse broadening than the energy-dependent broadening embedded within it. Secondly, 
we find that the correlations between E p and 5s are still positive and the correlations get 
much weaker after removing those merged pulses with high E p values. Therefore, if the 
merged pulses can cause high E p and the positive correlations between E p and 5s indeed 
exist the overlapping pulses may also make the 5s small. Based on the above analysis we 
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think the merged pulses would make the 5 values decrease. In order to identify further the 
correlations between E p and Ss more wider energy band and cleaner singe pulse data are 
needed. 

We further investigate the implications of the statistical relations. We find that the 
energy dependence of decay width (also the pulse width since it is dominated by the 
decay section) can be caused by the curvature effect+Band spectrum. Usually we study 
the observed light curve for a given energy or energy internal (e.g., BATSE has 4 energy 
channels). The light curve decays due to the delay arrival of high-latitude photons with 
lower flux, i.e., curvature effect. For a given observed energy less than the peak energy of 
Band spectrum, the intrinsic photon energy contributing to the energy will move from low 
energy part of Band spectrum to high energy one due to the lower Doppler factor at high 
latitude. This will naturally lead to the energy dependence of decay width and thus the 
pulse width. This in turn supports the decay phase is indeed due to the curvature effect, 
or at least related to it. In addition, different Band function parameters may lead to the 
different degree of pulse broadening. The main effect seems to be from the low-energy slope 
as analyzed above. The strong correlation between 5d and a also appears to result from the 
curvature effect and Band spectrum (see Fig 19]). 

The rise phase of the GRB pulse contributes less to the pulse width than the decay 
phase. However it is important and is considered to be relevant to the hydrodynamic time 
of shock crossing shell. Statistical studies have found that the rise phase is also related to 
the time-resolved Band spectrum in single pulse. For instance, some authors found that 
the time-resolved peak energy of pulse is tracking the pulse profile within a pulse for some 
bursts (Kaneko et al. 2006; Peng et al. 2009b, however, see Hakkila & Preece 2011 for 
different view). In current investigations, we didn't find correlations between broadening 
of the rise phase and the (time-integrated) Band spectrum parameters for different pulses, 
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suggesting that the formation of the rise phase is related to emission mechanism, while its 
energy dependence is independent of emission mechanism. The energy dependence of rise 
phase may arise from the intrinsic spectral evolution. As pointed by Fenimore & Sumner 
(1997), the intrinsic spectral evolution affects the rise phase much more significant than the 
decay phase. 

A deduction of the correlation between 5d and a we found is that the energy dependence 
of pulse with single power law (corresponding to a = (3 in Band spectrum) spectrum is 
systematically weaker than that with Band spectrum. If it is verified, this will support that 
the intrinsic Band spectrum and the curvature effect is indeed an important factor leading 
to the photon energy dependence of GRB pulses. Also it in turn supports that the decay 
phase of pulse is produced by the curvature effect. We will consider it in a future paper 
with Fermi/GBM or Konus-Wind data. 
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